State selection in the noisy stabilized Kuramoto Sivashinsky equation 
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In this work, we study the ID stabiHzed Kuramoto Sivashinsky equation with additive uncorre- 
lated stochastic noise. The Eckhaus stable band of the deterministic equation collapses to a narrow 
region near the center of the band. This is consistent with the behavior of the phase diffusion 
constants of these states. Some connections to the phenomenon of state selection in driven out of 
equilibrium systems are made. 
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I. INTRODUCTION 



In this paper, the old question of pattern selection in 
extended systems is studied in a simple but non triv- 
ial driven out of equilibrium model with additive Gaus- 
sian distributed stochastic noise. A number of physical 
systems fall into this category, examples of which are 
wavelength selection in directional solidification [1] and 
Rayleigh-Benard convection [2]. The absence of selec- 
tion is claimed in models of convection [3] and in direc- 
tional solidification [4, 5]. Contrary claims have also been 
made [6~9]. The apparently simple model we choose to 
study is the stabilized Kuramoto Sivashinsky equation in 
one spatial dimension with additive uncorrelated Gaus- 
sian distributed stochastic noise. This equation has the 
essential ingredient of non linearity, has a band of sta- 
ble steady states in the absence of noise, displays most 
of the instabilities and stationary states of real systems 
[10] and is sufficiently simple to allow a detailed analy- 
sis in the presence of stochastic noise. All real systems 
are subject to some sort of stochastic noise which may 
have very important consequences. For example, a sys- 
tem evolving towards equilibrium will ultimately end up 
in a unique stationary equilibrium state which is the state 
of minimum free energy. There may be several states cor- 
responding to local free energy minima but the system 
eventually reaches the state of absolute minimum free 
energy. Should the system evolve in a noiseless deter- 
ministic fashion, it will become stuck in the first local 
minimum encountered which implies that its final state 
depends on its initial state. However, when stochastic 
noise is present, the system will escape from any local 
minimum and will eventually end up in the absolute free 
energy minimum. This picture is consistent with selec- 
tion of the unique equilibrium state at very late times as 
in [11, 12]. 



The temporal evolution of many driven out of equilib- 
rium systems cannot be described in terms of a potential 
as 



dt 



+ c 



(1) 
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where C(r, t) is a Gaussian distributed uncorrelated ran- 
dom noise with (C(r,i)) = and (C(r, <)C(0, 0)) = 
2e6{r)d{t). A system whose evolution is described by 
Eq. (1) will reach the stationary state corresponding to 
the absolute minimum of the potential as t ^ oo, 

as in the Swift-Hohenbcrg equation [13]. The stochastic 
noise C(r,i) in Eq. (1) is essential for the selection of a 
unique stationary state provided that the noise strength 
e is not too large [12]. If we interpret the potential T{(j>) 
as the free energy, then e oc T, the temperature of the 
system. As expected, when e > Cc, or T > Tc, no unique 
state is selected and the system is disordered [12]. Also, 
selection of a unique equilibrium state occurs only in the 
thermodynamic limit when the probability of the state 
is unity. In a finite system, it is well known that fluctu- 
ations or noise ensure that all possible states are visited 
with finite probability so there is no true selection. 

The role of noise is not well understood, and mostly ig- 
nored in non-equilibrium systems on the grounds that it 
is so small that its neglect is justified [2]. An important 
question that has been addressed both experimentally 
and theoretically, is the effect of noise on the stability of 
stationary states or patterns[l, 14-21]. Pattern selection 
is a widely investigated phenomenon. Systems where this 
has been observed include convective rolls in Rayleigh- 
Benard instabilities [22-24], propagation of Taylor vor- 
tices in unstable Couette- Taylor flow [25], electrohydro- 
dynamic convection in nematic liquid crystals [26, 27], 
directional solidiflcation [28-30] and flngering instabili- 
ties in Hele-Shaw type experiments [31, 32]. These stud- 
ies show that fluctuations seem to play an important role 
and that they should not be ignored. In this paper we 
further examine this phenomenon in an attempt to ob- 
tain a better understanding of the effects of noise. 

In dissipative systems, spatially periodic structures 
emerge as a result of a primary instability which develops 
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when the control parameter exceeds its critical value [33] . 
The periodic structures can themselves become unstable 
and develop secondary instabilities under certain pertur- 
bations of the primary pattern [2]. For example, the Eck- 
haus instability arises from the translational invariance 
of the periodic solution. As a result of this instability, the 
system undergoes a change in the wavelength of the pat- 
tern. We consider the stabilized Kuramoto Sivashinsky 
(SKS) equation [10] 



dtu = -mt - V^it - V^w + (Vm)^ 



(2) 



where u{x,t) is a scalar function in one space dimension 
describing an interface profile. The control parameter, 
a > 0, is a stabilizing mechanism. The SKS equation is 
one of the simplest equations describing dissipative sys- 
tems and is used to describe directional solidification [10] 
and the Burton-Cabrera-Frank model of terrace growth 
[34]. Extensive studies [10, 35] reveal that this mod- 
est equation displays very rich dynamics and has a vari- 
ety of secondary instabilities: (i) Eckhaus instability, (ii) 
period-halving of the cellular state, (iii) parity breaking, 
(iv) vacillating breathing and (v) oscillation with a spa- 
tial wavelength "irrationally" related to the basic one, as 
well as an abundance of tertiary instabilities. A linear 
stability analysis of Eq. (2) about the trivial solution 

gives the dispersion relation 



w = with Su - e^*+*«^ 
A = —a + ~ showing that a primary Turing bifur- 
cation occurs at a critical value ac = 0.25 and ^ A=. 

v2 

For a < Uc, the trivial solution is unstable and modes 
in the band 1/2 + i/l/4 - a > q^ > 1/2 - y/l/4-a 
grow exponentially. However, the nonlinearity in Eq. (2) 
mixes these modes resulting in the emergence of stable 
periodic structures [10, 33, 35]. The stability diagram of 
the SKS equation of Fig. (1) shows the neutral curve and 
the different regions of secondary bifurcations. 

In this work, we study the effect of noise on static cells 
in the stable region inside the Eckhaus stable band. The 
paper is organized as follows. In Section II, we present 
the numerical simulations of the noisy SKS equation and 
the results obtained. In Section III we discuss the lin- 
ear stability analysis of the deterministic noiseless SKS 
equation and the phase-diffusion coefficients of the steady 
states. We also present numerical computation of the 
phase-diffusion coefficients of the steady states of the de- 
terministic SKS using Auto [36, 37], a software for con- 
tinuation and bifurcations of a system of ordinary differ- 
ential equations (ODE). We compare the results obtained 
from direct simulations of the noisy SKS equation with 
the computation of the phase-diffusion coefficient. Fi- 
nally, in Section IV we draw some conclusions from our 
study. 



II. THE NOISY KURAMOTO-SIVASHINSKY 
EQUATION-NUMERICAL SIMULATION 

The Eckhaus stable band of stationary periodic states 
for the deterministic SKS equation is shown in Fig. (1) 
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FIG. 1: Stability diagram of the primary and secondary states 
in the SKS equation. The domains of static cells, oscillating 
cells (OSC), drifting cells (DRIFT), and oscillating-drifting 
cells (OSCD), are bounded respectively by open circles, dark 
circles, open squares, and dark squares. In empty domains, a 
given state is unstable. The full line is the neutral curve of 
the mode of wave number q, and the dashed line for the mode 
2q. The figure is taken from Ref. [35] (with permission). 



[10, 33, 35]. All states inside this band are stable against 
small perturbations while those outside the band are un- 
stable against long wavelength perturbations and decay 
to a state inside the band which depends on the exact 
form of the perturbation. Since any real system is sub- 
ject to a variety of external and internal perturbations 
such as thermal ffuctuations, vibrations from a slammed 
door, a heavy truck outside the laboratory, etc, it is of 
some interest to study the consequences of random noise 
on these stable stationary states. The simplest way to 
represent such effects is to add uncorrelated Gaussian 
distributed noise to the deterministic evolution equation. 
In this section, we investigate numerically the effects of 
additive Gaussian distributed white noise on the cellu- 
lar stationary states of the SKS equation. The Langevin 
equation associated with Eq. (2) is 



dfU = 



-au-V^u-V^u+{\7u)^ + C{x,t). (3) 



where (^{x, t) is a Gaussian distributed noise with 
(C(a;,<)) = 0, and {C{x,t)C,{x' ,t')) = 2eS{x - x')S{t - t') 
with denoting an average over the noise distribution. 

The noisy SKS equation of Eq. (3) has been studied in 
the context of rough growth and morphological instabili- 
ties in various growth processes such as electrodeposition 
[40-42] and of the evolution of a surface undergoing ion 
sputtering [43, 44]. The focus of these earlier works is to 
understand the forced growth into a rough interface with 
emphasis on the interface width and its scaling proper- 
ties. In this work, we focus the formation of a periodic 
stationary state and the the role of the stochastic noise 
in the selection of a unique stationary state rather than 
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the scaling of the mesoscopic width of the interface. 

To perform numerical simulations, we discretize Eq. 
(3) by defining uj — uixi^tj) with Xi — iAx, tj = jAt. 
The spatial grid of N points is labelled by the subscript 
1 < i < N and the time grid of J\f points by the super- 
script I < j < Af. We impose periodic spatial bound- 
ary conditions u-'-^j^ = and. since we have an initial 
value problem and a first order time derivative, specify- 
ing u{x,t = 0) completes the definition of the problem. 
The discretized Ito-Langevin equation corresponding to 
Eq. (3) is [45] 



7-1-1 3 

U- = U; 



F/ \u\ = —aui - 



1 



2eAi 



(Ax)- 



-t- 



1 



4(Aa;)- 



6ul ~^u\_^+u\_^ 
uUf (4) 



where r\\ is a Gaussian white noise of unit variance with 
= ^iid (rifrjl.) = SikSji. Because of the term in Eq. 
(4), (y^u)l = {Ax)-^{ui^,-^^, + 6i4~4w>_, + w>_^ 
for 1 < i < TV, it is convenient to define a grid of iV -I- 4 
points labelled by a = -1,0,1 ■ ■ ■ N, N + 1, N + 2 and 

We restrict our study to stationary periodic states 
by limiting the control parameter to 0.16 < a < 0.25 
[10, 35]. In the discrete system, the allowed wavenum- 
bers are q = 2ttNc/Lx where 1 < A'^c < is an integer 
corresponding to the number of cells in the pattern and 
Lx = NAx is the size of the system. Most of our simula- 
tions are done with = 1024 and grid spacing Ax = 0.50 
which is the largest value for which we could discern no 
difference in the pattern by changing Ax Ax/2. When 
we compare results for Ax = 1 and Ax = 0.50, the dif- 
ferences are visually very obvious but not for Ax 0.50 
and Ax = 0.25. The error due to the spatial discretiza- 
tion is 0{Ax^). Because of the V^u term in Eq. (3), 
the Courant-Friedrichs-Lewy (CFL) condition determin- 
ing the convergence of the discretized form in Eq. (4) 
to the solution of the continuous PDE of Eq. (3) is 
At < C{Axy. Numerically, we find C « 0.1331 and 
we choose At = 0.006 which satisfies the CFL condition 
when Ax = 0.50. For Ax = 0.25 the CFL condition 
requires a time step At = 0(10"^) which implies an or- 
der of magnitude greater computational time Afr than 
for Ax = 0.50 for the same real time T = AfAt. Here, 
T oc A'' is the time to perform a single update t ^ t + At. 

Direct simulations on systems of finite size, = 
A^Ax, show that noise destabilizes all steady states of 
the deterministic SKS equation of Eq. (2), as expected. 
A truly stable stationary state can exist only in the ther- 
modynamic limit Lx ^ oo which implies that a finite 
size scaling analysis, similar to that used to study equi- 
librium phase transitions [46] , might provide a more con- 
vincing case for state selection but this is beyond our abil- 



ity. Thus, our conclusions from these simulations must 
rely on some operational definition of stability which we 
choose as the existence of the state in the presence of 
noise for at least some specified number of time steps. 
This measure allows us to compare the stability of differ- 
ent states. 

The Eckhaus band of spatially periodic stable states 
shrinks in the presence of noise. Even though we expect 
the stable band to shrink to a point, in the simulations we 
could only observe it to shrink to a narrow band in most 
cases. We were not able to further narrow down this 
region to a single state, because of limitations in both 
the computational time and allowed noise strength. In 
the absence of a clear definition of an appropriate noise 
strength, we limited ourselves to noise strengths for which 
we could still clearly observe a periodic pattern. The fur- 
ther from the boundaries of the Eckhaus band the funda- 
mental wavenumber q is chosen, the more stable against 
noise the periodic state becomes, so that a larger noise 
strength e, or, equivalently, a longer simulation time is 
required to destroy the state. As a is decreased, states 
inside the Eckhaus stable band become more stable. 

In the simulations, we take Ax ~ 0.50, At = 0.006 and 
A^ = 1024 since these values allow for the largest real time 
T = MAt while satisfying the CFL condition. A few sim- 
ulations with A^ = 2024 were done with no change in the 
results. For a — 0.24, the band of stable wavenumbers 
is 53 < A'c < 61 corresponding to 0.65 < Qs < 0.74. In 
the presence of noise with e = 1.10^^, the initial deter- 
ministic static steady state with Nc = 54 evolves into the 
periodic state with Nc = 55 shown in Fig. (2). Similarly, 
the Nc = 60 deterministic steady state transitions to the 
Nc = 59 periodic state. How long a state of a particular 
periodicity survives before transitioning to another peri- 
odicity depends on the noise sequence generated. Also, 
the lower the noise strength e, the longer a particular 
periodic state survives before transitioning to a neigh- 
boring periodicity. For example, in a typical simulation 
with e — 2.10""*, the system transitioned from Nc = 59 
to Nc = 58 after 7.10^ time steps. When e = 5.10""*, 
the same A^c = 59 state evolved first into an Nc = 58 
state and then to a A^c = 57 state in 3.10^ time steps. 
The Nc = 55 state evolves into an Nc = 56 state. The 
Nc = 56 and Nc = 58 states in turn evolve into an 
Nc = 57 or Qs ~ 0.6995 state which appears to be the 
most stable state as it remains stable for 10^ time steps 
and noise strength up to e = 5.10"^. 

As a is decreased, the states in the Eckhaus band be- 
come increasingly stable against noise and we could not 
shrink the Eckhaus band to a point but only to a nar- 
rower band. For a — 0.2, the Eckhaus band consists of 
0.58 < qs < 0.77, corresponding to 48 < A'c < 63. In 
the presence of noise, states in the middle of the band 
survive. In Fig. (3) is shown the evolution of an initial 
Nc = 61 state to an Nc = 59 state at lO*" time steps with 
noise strength e = 5.10"**. This state is not the most 
stable as six states with 53 < A'c < 58, corresponding 
to 0.65 < Qs < 0.71, remain stable for 1.10^ time steps 
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FIG. 2: Time evolution of an A'^c = 54 state at a = 0.24 in the presence of noise with e — 1.10 




100 200 300 400 500 600 
lattice position i 




150 
100 
50 


- ( 

150 
100 
50 


- ( 

150 
100 
50 


- ( 

150 
100 
50 







20 40 59 80 100 120 






20 40 59 80 100 120 






20 40 60 80 100 120 







100 200 300 400 500 600 
lattice position i 



FIG. 3: Time evolution of an Nc = 61 state at q = 0.2 in the presence of noise with e = 5.10 



against noise of strength e < 2.10 ^. 

For a ~ 0.17, Fig. (4) shows the evolution of a periodie 
state with Nc — 63 into a. Nc = 62 periodic state after 
1.10^ time steps with noise strength e = LID""*. Fig. (5) 
shows the evolution of a A'^c = 48 to a iVc = 49 periodic 
state with a larger noise of strength e = 5.10""*. In this 
figure, small amplitudes of periodicities Nc = 48 and 50 
with a main peak at Nc = 49 can be seen in the snapshot 
of the pattern at 7.10^ time steps which have become 
much smaller by 10^ time steps when this simulation was 
terminated. In Fig. (5), there is also a small second 
harmonic peak at 2Nc in the initial deterministic steady 
state with Nc 48 and in the Nc = 49 state after 10'^ 



time steps. The Eckhaus band at a = 0.17 contains 
periodicities 45 < iVc < 65 or 0.55 < Qs < 0.80 collapses 
into a narrower band of periodicities 53 < Nc < 57 or 
0.65 < Qs < 0.70 in the presence of noise. All states in 
this narrow band remain stable against noise of strength 
e = 2.10"^ for 1.10^ time steps. For a = 0.17, states 
near the center of the Eckhaus band are very stable and 
destroying them by direct simulation of Eq. (3) with 
noise strengths for which a periodic state clearly exists 
can take a very long computational time. 
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FIG. 4: Time evolution of an Nc = 63 state at a = 0.17 in the presence of noise with e — 1.10 
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FIG. 5: Time evolution of an A'^c ~ 48 state at a = 0.17 in the presence of noise with e = 5.10 



III. PHASE DIFFUSION COEFFICIENT 

A. Stability Analysis 

To understand the destabilization of the static steady 
states in the presence of noise, we linearize the profile 
u{x, t) about a periodic stationary state of period L, 
u*(x) = u*{x + L) 

u{x,t) ~ u*{x)+v{x,t) 

■ = -{a + d^ + d^K + {d^u*f 
vt = [~a~dl-di + 2id^u*)d^]v = C,v (5) 



where v is an arbitrary function v{x,t) G C", and 
is a linear operator which depends on the steady state 
solution u*{x). To check for the stability of the steady 
state, we compute the spectrum of >C*, which reduces to 
solving the eigenvalue problem 



C^v - Xv = 
We define a Bloch wave operator 

:= -a - {d, + - (a, + vf 
where G C. Eq. (6) becomes 

[a + {d., + vf + (a, + vf - 2(a,u*)(a, + 1^) + \\v = o(8) 



(6) 

2{d^u*)[d^ + v) (7) 
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FIG. 7: A|| = 2D{q) as a function of the period L for a — 0.2. 
The wavenumber q = 2n/L. 



for some v{x) = v{x + L) G C" and some e i[0, ^). 
We focus our attention on the eigenvalue Ao defined such 
that 5RAo > KAj Vj, which determines the stabihty of the 
mode u* {x). 

The Eckhaus instabihty is a long wavelength instability 
implying that = is the important spatial eigenvalue 
[2]. Expanding Ao(i') about = 0, we obtain 

Ao(i') - Ao(0) + ^^^l.=o + Y^I-=o + • • • (9) 

We define the group velocity Cg := — ^|^=o G 1^ and 

the phase-diffusion coefficient D := ^-d^\iy=a G IR- For 
stable static cells, Ao(0) = and Cg ~ 0. Thus, D{q) 
determines the stability of the periodic steady state with 
wavenumber q [38, 39]. 



B. Computing the Phase Diffusion Coefficient 

To find the eigenvalues as functions of v and the dif- 
fusion coefficient D{q), we use Auto [36, 37] which is a 
software for continuation and bifurcation of a system of 
autonomous ordinary differential equations (ODE). To 
use Auto, we rewrite the steady state equation associ- 
ated with Eq. (2) as a system of first order ODEs and, 
to get the correct number of parameters to solve the SKS 
equation, we add the term cu^, where the parameter c 
is the speed of the wave. This extra parameter is also 
used to check the results of the computations since it 
must always come out to be zero in our case. We write 
U = {u,Ux,Uxx,Uxxx) which allows us to normalize the 
spatial period L to unity. We thus consider the boundary- 
value problem on the interval (0,1) [36]. As a first step, 
we must find the stationary solutions of the deterministic 
SKS equation within Auto as we use these solutions as 
input to the stability calculations. The stationary solu- 
tions U obey the vector form of Eq. (5) 

U^ = LF{U,c)U (10) 

with the boundary and normalization conditions [36] 

U{1) = [/(O) 

/ {U'{x),Uoid{x)-U{x))dx^Q, 
Jo 

where Uoid(x) is the vector U from the previous step in 
the iterative computational procedure, U' = Ux and 

/ 1 \ 

mc)= I I J ; . (11) 

\-a (ux-c) -10/ 

Linearizing Eq. (10) about a stationary solution U* by 
writing U = U* +V with V = {v, Vx,Vxx,Vxxx) gives the 
vector form of Eq. (8) 

Vx = L[A + XB-:y]V 
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FIG. 9: qmax, the wave number maximizing D{q) as a func- 
tion of a. Error bars are largest possible errors. 



expected for stationary periodic patterns. The non van- 
ishing part of the output is A|| (L) as a function of the spa- 
tial period L where A|| = 2D{q), the phase diffusion con- 
stant for fundamental wave number q = in/L. D{q) < 
for periodic states outside the Eckhaus band indicating 
that these states are unstable against long wavelength 
perturbations. Following the standard interpretation, we 
identify the boundary of the Eckhaus band by the solu- 
tion of D{a,q) = and the stable band as the region 
D{a,q) > [2, 10]. From Figs. (6, 7, 8), we see that 
the maximum of A|| rises dramatically as a is decreased 
from 0.24 to 0.17. The wave number qmax maximizing 
X\\{q) = 2D{a, q) is shown as a function of a in Fig. (9). 
The error bars in Fig. (9) are an overestimate as they 
represent the maximum possible uncertainty in an ex- 
trapolation from nearby points from Auto. We interpret 
these values of qmaxioi) as the selected wave numbers as 
discussed in Section IV. 



A 



B = 







V-1 



1 











c) 




(12) 



We supplement Eq. (12) with the boundary and integral 
conditions 



V{Q) = V{1) 
j^{Void{x).V{x))dx = l. 



(13) 



Define Ai 



^1 ■= "dlTli'^O' ^^"^ '^ll ^d^\''=o ^nd denote T4 
by V' . Differentiating Eq. (12) with respect to v and 
evaluating the result for the most dangerous eigenvalues 
A = = gives f36l 



L[AV\ + (A|B - l)V] 

L[AV\\+2(\B -l)V\+\\BV] (14) 



We impose boundary and integral conditions 



{V{x),V\{x))dx = 0= / {V{x),V\\{x))dx 



M(o) = 
Mi(o) 



^(1) 



(15) 



and obtain A| and A|| by solving Eqs. (10) to (15) using 
Auto. To guarantee convergence, we use a finite differ- 
ence scheme to generate a steady state solution of the 
SKS equation, Eq. (10), and use this as input for solving 
Eqs. (12-15). The output of Auto consists of the eigen- 
values A|(L) and A||(L) as functions of the period L and 
the value of the wave speed c. A| and c are always zero as 



IV. CONCLUSIONS 

The stability of the stationary periodic states of the 
deterministic SKS equation is studied (i) by simulation 
of the SKS equation with additive uncorrelated Gaussian 
distributed noise and (ii) by computing the decay rates of 
perturbations about these periodic states. Direct simu- 
lations of the noisy SKS equation show that the Eckhaus 
boundary separating unstable and stable spatially peri- 
odic states of the deterministic SKS equation collapses 
in the presence of additive Gaussian distributed stochas- 
tic noise and the Eckhaus stable band shrinks to a much 
narrower band. For a Uc = 1/4, the simulations show 
that the Eckhaus band shrinks to a point as one of the al- 
lowed values of q is more stable than the others. However, 
as ttc — increases, the periodic states in the Eckhaus 
stable band become more stable against noise and the 
simulations could verify only that the width of the stable 
band shrinks. These results agree with the picture ob- 
tained from the deterministic equation by computing the 
phase diffusion coefficient D{q) of these spatially periodic 
states using Auto [36, 37]. 

We find that the value of the phase diffusion coefficient 
D{q) is negative for periodic states outside the Eckhaus 
band, positive for those inside and zero for the states 
at the boundary. This indicates that states inside the 
Eckhaus band are stable against long wavelength per- 
turbation. The phase diffusion coefficient seems to also 
determine the stability of the state against random fluc- 
tuations. The states with a higher value of the phase dif- 
fusion coefficient are more stable against noise. The max- 
imum of the phase diffusion coefficient curve increases as 
a decreases. 

For a = 0.24, simulations of the noisy SKS equations 
show that the most stable state in the presence of noise 
is Nc = 57 corresponding to g = 0.6995 with parameters 
N = 1024, Ax = 0.50. This state has the largest value of 
D{q) of the states with b'i < Nc< 61 comprising the Eck- 
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haus band for the N = 1024. Ax = 0.50, system used in 
our simulations. As a is decreased, the states inside the 
stable band become more stable against noise and, due 
to limitations in computational time and allowed noise 
strengths, we are able to narrow down this region to a 
single state for a ~ 0.24 only. However, the picture is 
consistent with that given by the phase diffusion coeffi- 
cient despite the large Aq = 2'k/{NAx) = 0.01227 ■•• be- 
tween allowed values of q because of the periodic bound- 
ary conditions. The stable region in the presence of noise 
for each a corresponds to a region located about the max- 
imum of the phase diffusion coefficient curve for that par- 
ticular value of a. Also, from Figs. (6, 7, 8), we see that 
the maximum value of the phase diffusion coefficient in- 
creases from about D « 2.75 at a = 0.24 to D « 9.5 at 
a = 0.17, which is consistent with the observation from 
direct numerical simulations that the stability against 
noise of states inside the Eckliaus band increases as a 
decreases. 

We note from our results from the phase diffusion 
coefficient of Fig. (9) that, as a ranges over the full 
range of stationary periodic states from a = 0.25 to 



a = 0.17, the periodicities vary over a very small range 
0.707 > qs > 0.685 corresponding to 58 > iV^ > 56 
when N = 1024 and Ax = 0.50. To obtain a larger 
range of Nc, larger system sizes are necessary. Perhaps 
N = 2^^ = 4096 is suitable but sizes like this will re- 
quire extremely long runs and a large increase in the 
noise strength e. In turn, this will make identification 
of the periodicity of a state problematic. Taken together 
these extra difficulties seem to preclude numerical testing 
of our hypothesis with our presently available computa- 
tional resources. 
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